New Tetrahedral Global Minimum for the 98-atom Lennard-Jones Cluster 
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A new atomic cluster structure corresponding to the global minimum of the 98-atom Lennard- 
Jones cluster has been found using a variant of the basin-hopping global optimization algorithm. 
The new structure has an unusual tetrahedral symmetry with an energy of — 543.665361e, which is 
0.022404e lower than the previous putative global minimum. The new LJgg structure is of particular 
interest because its tetrahedral symmetry establishes it as one of only three types of exceptions to 
the general pattern of icosahedral structural motifs for optimal LJ microclusters. Similar to the 
other exceptions the global minimum is difficult to find because it is at the bottom of a narrow 
funnel which only becomes thermodynamically most stable at low temperature. 
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The determination of the global minima of Lennard- 
Jones (LJ) clusters by numerical global optimization 
techniques has been been intensely studied in the size 
range A^= 13-147 by both chemical physicists and applied 
mathematicians |l]J^. The LJ potential, which is given 
by 
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(1) 



where e is the pair well depth and is the equilib- 

rium pair separation, is a simple yet reasonably accu- 
rate model of the interactions between heavy rare gas 
atoms. In general, there has been good agreement be- 
tween physical measurements on rare gas clusters from 
electron diffractometry [|j and mass spectrometry [||J|] 
and computational global optimization results regarding 
magic number sizes and corresponding cluster geometries 

. Both approaches find that Mackay icosahedra are 
the dominant structural motif. 

The L J microcluster problem has also become a bench- 
mark for evaluating global optimization algorithms. The 
number of local minima (excluding permutational iso- 
mers) on the potential energy surface (PES) is believed to 
grow exponentially with N [g,^ and is estimated to be of 
the order of 10'"' for A^=98. A wide variety of global op- 
timization techniques including simulated annealing , 
genetic algorithms |iri| - |l3|, smoothing and hypersurface 
deformation techniques |ll4P^ i lattice method s P ,P^, 
growth sequence analysis |l7|]l8|] , and tunneling |19| have 
been applied to the problem. Unbiased methods that 
make no assumptions regarding cluster geometry are of 
the most interest, since these have the best chance of suc- 
cessful generalization to more complex potentials such as 
those in the protein folding problem. 

Most of the global minima in this size range were first 
found by Northby in a lattice-based search of icosahedral 
structures m . These structures consist of a core Mackay 



icosahedron (Figure |l|b) surrounded by a partially filled 
outer shell. More recently, there have been a number 
of improvements in some of these putative global min- 
ima. Firstly, further refinements to Northby's algorithm, 
particularly the relaxation of the assumption that the 
core Mackay icosahedron is always complete, has a led 
to a number of new global minima [ p^yi^ , |20| , ^ . Sec- 
ondly, consideration of particularly stable face-centred- 
cubic (fee) and decahedral forms has also led to new 
global minima ||l9|,^,^ . At A^=38 the global minimum 
is a fee truncated octahedron (Figure |^a) and at N—75- 
77 and 102-104 the global minima are based on Marks 
decahedra (Figure |l|c) . Thirdly, powerful unbiased global 
optimization algorithms, particularly the basin-hopping 
and genetic algorithms |12 1^, have recently begun 
to catch up with those methods that incorporate partic- 
ular physical insights into the LJ problem, and are now 
able to find all the known lowest-energy minima. 




FIG. 1. Three particularly stable examples of the known 
morphologies for LJ clusters: (a) the 38-atom fee truncated 
octahedron, (b) the 55-atom Mackay icosahedron, and (c) 
the 75-atom Marks decahedron. 

Given this combined attack on the LJ optimization 
problem, it might have been imagined that all the global 
minima for N < 150 had been found. Here, however, 
we report a new lowest-energy structure for LJgg, which 
has an energy of — 543.665361e and Tfi point group sym- 
metry. This compares to an energy of — 543.642957e for 
the previous icosahedral putative global minimum which 
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was found by Deaven et al The LJgg global mini- 

mum is organized around a central foe tetrahedron with 
four atoms on each edge (Figure ||c). Four additional fcc 
tetrahedrons (minus apices) are erected over the faces 
of the central tetrahedron to form a 56-atom stellated 
tetrahedron (Figure ||b). An additional 42 atoms dec- 
orate the closed-packed sites on the surface of the stel- 
lated tetrahedron to complete the structure (Figure |2^) . 
The new LJgg structure is of particular interest because 
its tetrahedral symmetry establishes it as only the third 
known type of exception to the general pattern of icosahe- 
dral structural motifs for optimal LJ microclusters, and 
the first to be discovered by an unbiased optimization 
method. 

Given its unusual structure one might wonder why it 
is so low in energy. For LJ clusters optimizing the energy 
is a balance between maximizing the number of nearest 
neighbours and minimizing the strain energy (the ener- 
getic penalty for nearest-neighbour distances deviating 
from the equilibrium pair value) . The spherical shape 
and high proportion of {111} faces gives the structure a 
large number of nearest neighbours (432 compared to 437 
for the lowest-energy icosahedral minimum and 428 for 
the lowest-energy decahedral structure), whilst its strain 
energy is intermediate between icosahedral and decahe- 
dral structures. The lower strain energy allows it to be 
lower in energy than the icosahedral minima, even though 
it has fewer nearest neighbours. The strain in the struc- 
ture is focussed around the six edges of the central fcc 
tetrahedron. The atoms along these edges have the same 
local coordination as atoms along the five-fold axis of a 
decahedron. 




FIG. 2. (a) Front and back views of the new LJgg global 
minimum, (b) The 56-atom stellated tetrahedron and (c) the 
20-atom tetrahedron that are at the centre of this structure. 



It is also natural to ask how general this structure is. 
Firstly, analogous structures can be formed with smaller 
and larger tetrahedra at their core. The previous one in 
this series is at iV=34 and the next one is at 7V=195. 
However, these structures are not energetically compet- 
itive: the former because it has too high a proportion 
of {100} faces, and the latter because it is not sufh- 
ciently spherical. Secondly, the structures of the other 
non-icosahedral LJ global minima have been experimen- 
tally observed for gold and nickel clusters, and 
found to be particularly stable in theoretical calculations 
of transition metal clusters |^^. Therefore, we performed 
some optimization calculations for the Sutton-Chen fam- 
ily of potentials |^ . The tetrahedral structure was low- 
est in energy for silver, but a decahedral minimum was 
lower in energy for nickel and a fcc minimum for gold. 
This is consistent with previous results, which indicted 
that, of these three metals, silver clusters exhibited or- 
dered structures with the most strain ||27|| . 

The new LJgg optimum was found using a variant 
of the basin- hopping global optimization algorithm pi] . 
The key idea behind the algorithm is the mapping of 
the original LJ potential energy function, £'(x), for each 
point X on the 3A'^-dimensional Cartesian coordinate 
space onto a "transformed" energy function, r(x). T(x) 
takes the value of -^(x) at the local minium, Xmin, ar- 
rived at by applying a given local optimization proce- 
dure, such as the conjugate gradient algorithm, with x 
as the starting point for the algorithm. Thus T(x) is 
a "plateau" function that takes on the constant value 
-£'(xmin) on the catchment basin surrounding each local 
minimum Xmin- T'(x) is a lower bound to i?(x) and co- 
incides with -E(x) at all of the latter's local minima, but 
all barriers are removed in the T'(x) landscape and tran- 
sitions bewteen catchment basins can take place all along 
the basin boundaries. 

The original basin-hopping algorithm consists of a 
Metropolis search of the transformed landscape, r(x), 
using a Monte Carlo sampling procedure to move be- 
tween local minima. In the variant used in the discovery 
of LJgg 1^^, the Metropolis criterion of accepting up- 
hill moves with a probability that is an exponentially de- 
creasing function of the energy increment is abandoned 
in favor of only accepting downhill moves. The algorithm 
is restarted from a fresh random starting local minimum 
whenever progress stalls for a sufhciently large number 
of move attempts. The variant was successful in locat- 
ing the LJgg global minimum in 6 of 1000 random starts, 
with a mean computational time between encounters of 
about 30 hours on a 333 MHz Sun Ultra II processor. 
This structure has also been subsequently found using 
the original basin- hopping algorithm [^0| . 

These results show that the LJgg global minimum is 
particularly difficult to find. The origins of this diffi- 
culty are probably similar to the other non-icosahedral 
clusters. Analyses of the PESs of LJag and LJ75 using 
disconnectivity graphs have shown that they consist of a 
wide icosahedral "funnel" l3^,p3| and a much narrower 



2 



funnel leading to the global minimum . On relax- 

ation down the PES the cluster is much more likely to 
enter the icosahedral funnel, where it is then trapped be- 
cause of the large (free) energy barriers to escape from 
this funnel into the funnel of the global minimum. 

This situation is compounded by the thermodynam- 
ics of these clusters js^,^^. The icosahedral funnel has 
a larger entropy because of the larger number of low- 
energy minima, and so the funnel of the global mini- 
mum is only lowest in free energy at low temperatures. 
Therefore, at temperatures where the dynamics occur at 
a reasonable rate there is a thermodynamic driving force 
to enter the icosahedral funnel. For LJgg there are at 
least 114 minima that are lower in energy than the sec- 
ond lowest-energy minimum in the tetrahedral funnel, 
and so the global minimum is only lowest in free energy 
below T=0.0035eA:~^ (a typical melting temperature for 
a LJ cluster is 0.3ek~^). This transition t emp erature is 
markedly lower than for LJ38 |Q or LJ75 |2J. 

The basin-hopping transformation of the PES helps to 
ameliorate some of these difficulties. The transformation 
changes the thermodynamics so that the global minimum 
still has a significant occupation probability at tempera- 
tures where the cluster can escape from the icosahedral 
funnel. However, on relaxation down the PES the system 
is still much more likely to enter the icosahedral funnel. 
For example, our optimization runs were fifteen times 
more likely to terminate at the lowest-energy LJgg icosa- 
hedral minimum than at the global minimum. 

Coordinate files for the new LJgg structure, as well as 
all other putative LJ microcluster global optima, can be 
found in the Cambridge Cluster Database p7|. 
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